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ABSTRACT 



The field-induced dynamics of artificial spin ice are determined in part by interactions 
between magnetic islands, and the switching characteristics of each island. Disorder in either 
of these affects the response to applied fields. Numerical simulations are used to show that 
disorder effects are determined primarily by the strength of disorder relative to inter-island 
interactions, rather than by the type of disorder. Weak and strong disorder regimes exist 
and can be defined in a quantitative way. 
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Artificial spin ice*' 3 is a two-dimensional frustrated system, consisting of arrays of single- 
domain sub-micron islands of magnetic material, which are coupled by stray magnetic fields. 
It was conceived as an experimental model of three dimensional spin ices^ such as Dy 2 Ti 2 07 
and Ho2Ti 2 07, with the advantage that its microstates can be imaged directly due to the 
relatively large island volume. A consequence of this large volume is that artificial spin ice 
is typically athermal, and dynamics must be induced by an external magnetic field, which 
acts to "tilt" the complex energy landscape and allow transitions between configurational 
states. 

Disorder in the system affects field-driven dynamics via two mechanisms. Disorder in 
island interactions (energetic disorder) alters the energies of configurational states. On the 
other hand, disorder in island magnetization switching characteristics (switching disorder) 
changes the barriers between states, altering the pathways the system takes through its 
space of configurational states^. There are many possible sources of quenched disorder. 
In this paper we study four disorder types and show that their effects are similar. We also 
distinguish two regimes of disorder strength, namely the weak disorder regime and the strong 
disorder regime. 

Previous experimental studies of disorder have shown that the hysteresis of artificial spin 
ice reveals the importance of disorder 6 . The switching characteristics of connected^ and 
disconnected 9 kagome ice arrays have been used to estimate the strength of disorder in 
those systems. In simulation studies of an analog to artificial spin ice based on vortices in 
patterned superconductors, disorder is shown to nucleate grain boundaries between different 
ground state orderings 4 . 

The four disorder types studied here are disorder in island positions, disorder in island 
orientations, a random perturbation to pairwise interactions and a random perturbation to 
island switching fields. Pairwise energy disorder is a simplified 'representation' of energetic 
disorder - in particular, it does not lead to correlations in island pair energies, unlike po- 
sitional and orientational disorder. In all cases, disorder is assumed to have a Gaussian 
distribution. For example, positional disorder is implemented by perturbing the x and y 

(i) 

coordinates of each island i by 8x' y , where the 5 come from a Gaussian distribution with 
mean and standard deviation o. 

We limit our considerations to square artificial spin ice 1 . In square geometry, four islands 
meet at each vertex of the array, as shown in Fig. [TJ The interactions of the four islands are 
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not all equivalent: four of the six pairwise interactions are stronger than the other two, due 
to the difference in distances and orientations between islands. As a result, there are four 
energetically-degenerate types of vertex configuration, and the ground state of the system is 
a tiling of Type 1 vertices. The population fraction of Type 1 vertices, rti, is a good measure 
of the energy of the system. 

In the absence of disorder, the four vertex types have distinct, single- valued energies 
Ei < E 2 < E 3 < E4. However, when energetic disorder is present the energy of a vertex 
depends on the particular disorder realization. Considering many disorder realizations leads 
to a distribution of energies for each vertex type. The width of the distribution enables us 
to characterize the strength of disorder. 

We make this characterization in a model system in which the island magnetizations are 
Ising point dipoles, interacting so that the dipolar energy of island i is 



with the island magnetization M and the nearest-neighbor distance both set to unity, so that 
the nearest-neighbor coupling has strength 1.5, in units of l/(47r/io). In a perfect system, 
this gives E x = -6 + y/2, E 2 = -y/2, E 3 = 0, £ 4 = 6 + y/2. 

Shown in Fig. [2^a-c) are energy bands for each vertex type, defined by the mean en- 
ergy plus/minus one standard deviation, plotted against disorder strength a for the three 
energetic disorder types. The bands become wider as the disorder strength increases, and 
eventually overlap occurs. We take the weak disorder regime to be that where there is no 
overlap, and the threshold of the strong disorder regime being the point where two bands 
first cross (in practice, the first crossing is always of the Type 2 and Type 3 bands, since 
they are closest initially). The transition occurs at <y(E pSuiT ) w 0.225 for pairwise energy 
disorder, <y(4>) ~ 0.3 for island orientational disorder and c(rj) ~ 0.05 for island positional 
disorder. (Note that = r XiV refers to each component of an island's position.) 

The fourth type of disorder we discuss has a different nature, because it is disorder in 
island magnetization switching characteristics. We model the switching of the Ising moment 
of island i as occuring when 





(1) 




rhi < — h, 



(2) 
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where h[^ t is the sum of the external field h and the dipolar field h# and m ; is the direction 
of the island's magnetic moment. 

Unlike energetic disorder, switching disorder leaves the energies of configurational states 
unchanged and cannot be characterized in terms of vertex energies. However, an analogous 
measure can be constructed by considering the applied field required for two single spin flip 
processes. When all switching fields are equal, the applied field required to convert a Type 
2 vertex into a Type 3 vertex is h c + l/ \/2, while the field required for the reverse process is 
h c — l/\/2. These processes are shown as insets to Fig. Eld). We take the boundary between 
weak and strong disorder to be the meeting point of the plus/minus one standard deviation 
bands for the fields required for the two processes. The standard deviation of the switching 
fields, a(h c ), is equal to the standard deviation in the applied field required, so the crossing 
point can be easily determined to be a = 1/V2 w 0.7, as seen in Fig. Eld). Note that this 
measure of disorder strength is relative to the inter-island coupling, not the mean switching 
field, since changing h c does not affect the disorder strength at which the bands meet. Thus, 
both energetic and switching disorder are measured relative to the same quantity and can 
be compared meaningfully. 

To see that the strength of disorder is more important than its origin, we compare in 
Figures [3] and H] the final Type 1 population of a 20 x 20 island square ice, attained under 
a rotating applied field for both a perfect system and systems subject to different types of 
weak disorder (Fig. [3]) and strong disorder (Fig. HJ). In the simulations, the field rotates 
with constant amplitude and angular step d6 = 0.01 radians for 10 cycles, long enough to 
obtain a steady state. At each field application, the system relaxes by flipping single spins 
according to criterion fl2]) until no further flips are possible. 

We have previously shown— that, in the absence of disorder, dynamics that create Type 
1 vertices can only occur for a narrow range of field amplitudes Ah ~ 2 and dynamics 
always start at array edges. Smaller fields are unable to flip spins, and larger fields force the 
magnetization to track the field. Between these limits, there are two field regimes. In the 
low field regime, 10.25 < h < 11.25, dynamics proceed in a regular process of spin flips that 
"invade" from the array edges. In the high field regime, 11.25 < h < 12, more dynamical 
processes are possible, but the final Type 1 vertex population is limited because dynamics 
must start at the array edges, and Type 1 vertices that are created near the array edges 
"block" this from happening. 
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The four different types of disorder studied all have similar effects on the n\ vs h curves. 
In particular, switching disorder has the same effect as energetic disorder. This points to the 
fact that dynamics are determined through the spin flip criterion (j2J) , which is affected equally 
by changes to the energy of configurations (through h^ ip ) and changes to the switching fields 
Ki . We expect similar results to hold for other switching criteria, such as Stoner-Wohlfarth 
switching, since the island magnetization reversal always depends on both coupling to other 
spins and the intrinsic properties of each island. 

On the other hand, even though different disorder types have similar effects, the difference 
between weak and strong disorder is striking. As expected, in the weak disorder regime 
disorder acts as a small perturbation on the perfect system. The two non-trivial field regimes 
are still observed, and Ah is changed little from the perfect system. In contrast, in the strong 
disorder regime the n\ vs h curve is strongly altered from the perfect system. Ah increases 
to 4: at very low fields, spins with unfavourable hai P or small he can flip, whereas at very 
high fields the magnetization cannot completely track the rotating field due to pinning of 
spins with favourable h,a p or large h c % \ The other significant change is that n\ is reduced 
in the low field regime and increased in the high field regime. This is a result of dynamics 
being allowed to start in the bulk. In the low field regime, this blocks the regular invasion 
process that would otherwise lead to large n\. In the high field regime, this allows Type 
1-creating dynamics to continue longer, increasing n\ over the perfect case. 

In conclusion, different types of disorder in artificial spin ice have very similar effects on 
simulated dynamics and can all be compared to an energy scale set by the nearest-neighbor 
coupling. We show elswhere^ that this allows disorder in an experimental system to be 
measured in terms of an effective switching field disorder. 

Z.B. and R.L.S. acknowledge the Australian Research Council for funding. Z.B. acknowl- 
edges funding from INFN and the Hackett Foundation. 
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FIG. 1. The sixteen vertex configurations of square artificial spin ice can be classified into four 
types, based on energy. 



a Disorder in spin couplings b Disorder in orientations 




FIG. 2. The mean energy plus/minus one standard deviation of a (top-to-bottom) Type 4, 3, 2, 
1 vertex, subject to (a) pairwise energy, (b) orientational and (c) positional disorder. The point 
where the Type 2 and Type 3 bands cross is taken to be the transition between the weak and strong 
disorder regimes. The distributions are calculated over 200 independent disorder realisations, (d) 
The mean plus/minus one standard deviation of the external field required to convert a Type 2 to 
a Type 3 vertex (upper band) and the reverse (lower band). The insets show the two processes. 
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FIG. 3. Final Type 1 population vs applied field h for a perfect 20 x 20 system (black open 
circles) and four different types of weak disorder: island position disorder with standard deviation 
in each component of the positions of 0.01 (green filled triangles); island orientation disorder with 
standard deviation 0.06 (red filled circles); island pair energies perturbed with standard deviation 
0.045 (blue filled squares); and switching field disorder with standard deviation 0.14 (orange empty 
triangles). In all cases, the disorder strength a = 0.2<r*, where a* is the transition from weak to 
strong disorder. 
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FIG. 4. Final Type 1 population vs applied field h for a perfect 20 x 20 system with open edges 
(black open circles) and four different types of strong disorder: island position disorder with stan- 
dard deviation in each component of the positions of 0.06 (green filled triangles); island orientation 
disorder with standard deviation 0.36 (red filled circles); island pair energies perturbed with stan- 
dard deviation 0.27 (blue filled squares); and switching field disorder with standard deviation 0.84 
(orange empty triangles). In all cases, the disorder strength a = 1.2a*, where a* is the transition 
from weak to strong disorder. 



9 



